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Abstract. We regard the Casimir energy of the universe as the main contribution to the 
cosmological constant. Using 5 dimensional models of the universe, the flat model and the 
warped one, we calculate Casimir energy. Introducing the new regularization, called sphere 
lattice regularization, we solve the divergence problem. The regularization utilizes the closed- 
string configuration. We consider 4 different approaches: 1) restriction of the integral region 
(Randall-Schwartz), 2) method of 1) using the minimal area surfaces, 3) introducing the weight 
function, 4) generalized path-integral. We claim the 5 dimensional field theories are quantized 
properly and all divergences are renormalized. At present, it is explicitly demonstrated in the 
numerical way, not in the analytical way. The renormalization-group function (/3-function) is 
explicitly obtained. The renormalization-group flow of the cosmological constant is concretely 
obtained. 



1. Introduction 

At present, the energy budget of the universe is not explained theoretically. We do not 
understand the dark matter and the dark energy. The latter one, which is the present main 
topic, is regarded as the same problem as the cosmological constant one, which is the long- 
lasting problemfTJ [2]. As for the problem, Polyakov[31 Sj made a notable conjecture that the 
cosmological constant, just like the QED coupling, flows to a small value in the IR region 
(screening phenomena). He made another comment [5 J that the dark energy, like the black body 
radiation 150 years ago, hides secrets of fundamental physics. These views about the problem 
look to be revived in the recent strong trend of the AdS/CFT (holography) approach to the 
condensed matter physics or to the viscous fluid dynamics(Navier-Stokes equation) [7] . 

The difficulty of the problem is strongly related to the unsuccessful situation of the quantum 
gravity. It also has the long research history since Feynman initiated in 1963[8j until the 
present string or D-brane research. In the microscopic side, the gravitational interaction has, 
when quantized, serious UV-divergences[9], while , in the macroscopic one, it has serious IR- 
divergences around the horizon (boundaries) |10t H] . 

About 2 years ago (2010 January), Verlinde|ll| made a shocking view about the role of the 
gravitational force. He made a very thoughtful analysis, including some Gedanken-experiments, 
about the gravitational force and finally reached the conclusion that the gravitational 
interaction (force) is not fundamental but is emergent (by some still-obscure mechanism). He 
referred some (not few) past literature about similar stand points such as Jacobson's[12] and 
Padmanabhan's[13]. The analysis demands very delicate treatment of the horizon which is 



the boundary of the theory. The physical quantities (IR-)diverge at the horizon. It strongly 
indicates the regularization problem in the near-horizon treatment where Hawking radiation 
(thermalization) occurs. 

The biggest discrepancy, in all physical quantities, between observation and theory appears 
in the cosmological constant A. 



S = J d 4 X^{^^-(R - 2 A) + LmatteA 



A/ " ~ N D 2 L , N DL = M pl R cos ~ 6. x 10 BU , (1) 



^obs 

where M p i (= and R CO s (= ^/Hq,Hq : Hubble constant) are Planck mass and the size 

of the universe respectively. N^l is an analog of Dirac's large number [2]. The cause of 
the theoretical difficulty lies in the two extreme-ends of mysterious branches of physics: the 
quantum gravity and the cosmology. 

Casimir energy is generally the vacuum energy of the free-part of the system dynamics. 
For the harmonic oscillator, it is the energy of the zero-point oscillation. By definition, it is 
independent of the coupling. It depends, however, on the boundaries and the system topology. 
It is caused by the quantum fluctuation. Because we have to deal with serious UV and IR 
divergences, highly-delicate regularization is required. The well-known one is that for the 
4D electromagnetism (free wave theory). The system of two parallely-placed metalic plates 
separated by the length I, has Casimir energy as follows. 

4D Space-time : V Ca .(!) = |p *b«(0 = ~ = 3^, B 4 (4-th Bernoulli No) = 
5D Space-time : V Cas (l) = - ^ C [f , F Cas (l) = ~ = -g^^T, C(5) = 1.03693- • • . (2) 
2. Casimir Energy of 5 Dimensional Electromagnetism 

For the electromagnetism in the flat 5D space-time (ds 2 = n^dx^dx" + dy 2 , n^ = 
diag(— 1, 1, 1, 1)), Casimir energy is given by [15J 

d 4 p 



E Cas (A,l) = / <A 7^4 I dy(F;(p,y)+4F+(p,y)) 



lp<\ (27r) 

. f°° ~ =F cosh k(2y — I) + cosh kl 
Ff(p,y) = -I dk± : y , J 1S , (3) 



v 



2sinh(A;Z) 



where =F is the y-parity (y o — y), I is the periodicity y — > y + 2l, and A is the momentum cut-off. 
p is the magnitude of the 4D Euclidean momentum (pi,P2,P3,P4: = ipo)- This is consistent with 
Ref.[IS]. 

As for the warped (AdSs) geometry, ds 2 = (r\ iliV dx^dx v + dz 2 )/uj 2 z 2 (uj: 5D curvature), the 
corresponding one is given by [T7] 



/n r I 1 r t ~*~ > 

/ dz FT(p,z),FT(p,z) = / {GT( z ,z)}dk 2 , 

GT( f) _ 2 / 2 {Io(|)Ko(pz)TKo(g)I (^)}{Io(|)Ko(pzQ T K (|)Io(pV)} 
P 2 * Z Io(f)K (|)-K (f)Io(|) 

1 The original Dirac's definition is [electromagnetic force] /[gravitational force] — G J^ c m ~ 2.3 x 10 39 (a = 

e 2 /4neohc) where the force is that working between the electron and the proton in a H-atom. The substitute, 
in the present case, is naturally [size of the universe] /[Planck length]. Note that Dirac's value is approximately 
equal to ([Planck mass] /[mass of nucleon]) 2 . 



where Io and Ko are O-th modified Bessel functions. For simplicity, we take the bulk particle 
mass M in such a way that M 2 = —Auj 2 < 0. T is defined by T = u;exp(— loS). 



3. Behavior of Casimir Energy 

The integral region, for the flat case ([3]), is shown in FigfTJ [i and e are the IR cutoff of p-axis and 
the UV cutoff of y-axis respectively. Figj2]shows the integrand of ([3]): p 3 (F7 +AF^) = p 3 F(p, y). 
In FigJ3] the behavior of the integrand of F~ is shown. The (inverse) table shape says the 




" Rayleigh- Jeans" dominance because Casimir energy density is proportinal to the cubic power 
of p in the region p <C A. Numerically Ec a s(A, is obtained as 



2tt 2 

(2tt)' 



-0.1249/A 5 - (1.41,0.706,0.353) x 10~ 5 I A 5 ln(ZA) 



(5) 



In FigfTJ the region below the hyperbolic curve is that of Randall-Schwartz(RS) [IS] . They 
claimed the rectangular region should be replaced by this restricted one for the purpose of 
reducing the divergences. 
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dy q 3 F(q,y) 



2tt 2 



-[-8.93814 x 10~ 2 A 4 



(6) 



which shows slightly milder than (0). ln(/A)-term does not appear. 

As for the Warped model, the situation is similar. FigJ3]shows the integral region. The z-axis 
range is u> _1 < z < T . We take IR-regularization-point of p as fi = AT/lu. In FigJS] the 
behavior of the integrand of F~^ as Q is shown. Figj6] shows the warped version of Figj3j E^ as 
@ is numerically obtained as 



2tt 2 
(2tt)' 



-0.0250- 



A 5 



(7) 




which does not depend on u and has no ln(A/T) term. When restricted to RS-region, the above 
value changes to 



2tt 2 A 5 f 9 4l A1 

— — r— <^ -1.58 x 1(T 2 - 1.69 x 1(T 4 In - } , (8 



which is independent of T and has the ln(A/o;) term. Divergence behavior is the same as the 
unrestricted case (|7j). 



4. Sphere Lattice Regularization 

The Randall- Schwartz's way of restricting the integral region does not sufficiently work to reduce 
the divergences. In Ref.[l9], we proposed a new way of the restriction based on the isotropy 
property of the system and the minimal area principle. Let us introduce two 4D hyper-surfaces, 
Buv and B/#, in the 5D bulk space-time. 

EW : r 1 = vV) 2 + (* 2 ) 2 + (z 3 ) 2 + (* 4 ) 2 = ruv(v) , e = j<y<l , 

B IR : p- 1 = y/(J)* + (x2)2 + ( X 3)2 + ( x 4)2 = rjfl(y ) ; £ = I. < y < J , (9) 

where x a (a = 1,2,3,4) is Euclidian coordinate. See FigJT] and FigjHJ The surfaces, in the 
"brane" located at y, are 3D sphere with the radius r{y). The radius changes along the extra 
axis y. The function form of r{y) is given by the minimal area (of the hyper surface) principle. 
We integrate the region bounded by B//j from below and by Buy from above as shown in FigJT) 
Some regularization parameters (A,fj,,A',fi',e) are defined there. The renormalization group 
interpretation is shown in FigjHJ The surface B[/y is stereo-graphically shown in Figj9j It shows 
the present regularization utilizes the closed-string configuration. The integral region for the 
warped case is shown in FigllOl 




Figure 4. Space of (z,p) for the integration 
(HJ. The hyperbolic curve is used in (jSJ). 



Figure 5. Behaviour of (— l/2)p 3 T (p, z) in 
(JH). T = l,w = 10 4 ,A = 10 4 . 1.0001/w < 
z < 0.9999/T, AT/uj < p < A. 



Figure 6. Behavior of In | \ F~ (k, z) \ = \n\k G" k (z,z)/(u;z) 3 \. uj = 10 4 ,T = 1,A = 2 x 10 4 . 
1.0001/w < z < 0.9999/T. AT/u < k < A. Note In | (1/2) x (1/2) | rj -1.39. 




5. Introduction of the Weight Function W(p,y) or W(p,z) 

Another way to reduce the divergences is, instead of restriction, to introduce the weight function 
W(p, y) or W(p, z) in the integration over 5D space-time. The determination of the function 
form is explained later using the minimal area principle. Now we assume some forms as sample 
examples, and check whether Casimir energy is finitely obtained. For the flat case, we take the 
following forms. 



W(p, y) = 



r nv (y) 



Figure 8. Regularization Surface Bm and 
Figure 7. Space of (p,y) for the integration b uv in the 5D coordinate space (x»,y), Flow 
(present proposal as the substitute of FigHJ. f Coarse Graining (Renormalization) and 

Sphere Lattice Regularization. 

p -t 




Figure 9. UV regularization surface in 5D 

coordinate space. Figure 10. Space of (p,z) for the integration 

(present proposal as the substitute of FigH]). 



(iVO-ie^ 1 / 2 )^ 2 -^/^ 2 /' 2 = Wi{p,y), N x = 1.557/8^ 2 elliptic suppression 

(iY 2 )~V^ = W 2 (p,y), N 2 = 2(lA) 3 /8ir 2 hyperbolic suppression! 



Wi is considered to imitate the Randall-Schwartz restriction. FigJl II shows Wi-weighted case of 
Figj2j We notice, in the "valley", the depth, the location and the shape change. The Casimir 



Figure 11. Behaviour of p 3 W\(p, y)F(p, y) (elliptic suppression). A = 10, / = 1 
1/A < y < 0.99999/, l/l <p < A. 




energy is numerically obtained as 



^Cas 



f -(2.500, 2.501, 2.501)^, + (-0.142,1.09,1.13) x 10~ 4 Ab ^ A) for Wi , . 

\ -(6.0392,6.0394,6.03945) x 10- 2 ^ - (24.7,2.79,1.60) x 10- 8 ^|^ for W 2 

Triplet data show unstableness of the value due to insufficient range of I and A. ln(ZA) term 
vanishes for W 2 in the present calculational precision. 

For the warped case, the situation goes similarly. We take the following weights. 

W(p,z) = 

{ (iV^-V^PV^-a/^r 2 = W - 1 (p jZ ) > #1 = 1.711/87^ elliptic suppr. 

\ {N 2 )~ 1 e-P zT /" = W 2 (p,z), N 2 = 2^/8vr 2 hyperbolic suppr. 1 l ~ 

Figll2l is Wi-weighted case of Figj5j The "valley" changes in its depth, shape and location. 
Casimir energy is given by 

f *£Ax 1.2/l + O.H In £-0.10 In 4) for Wi 
Cas I ^ A 4 x 0.062 1 1 + 0.03 In £ - 0.08 In A} for W 2 

For the hyperbolic case (W2), the divergence is slightly milder than the case of RS-restriction 
© (A 5 -+ A 4 ). It is new that ln(A/T) appears. 

After calculating for 13 different weights [T5"l 1 17 j . we conclude Casimir energy can be expressed 
as follows (except the hyperbolic weights). For the flat case, 

£cL/AZ = -^(l-4cln(ZA)) . (14) 



Figure 12. Behavior of (— Ni/2)p 3 W\(p, z)F (p, z) (elliptic suppression). A = 20000, oj = 
5000, T = 1 . 1.0001/w < z < 0.9999/T, /z = AT/u < p < A. 




For the warped case 

E^/Ar 1 = -aw 4 (l-4cln(A/w)-4c'ln(A/r)) . (15) 

The parameters a, c, c! depend on the form of the chosen weight. The factors Al and AT -1 are 
the area of the rectangular regions (normalization constants). This result says the divergences of 
the 5D Casimir energy reduces to the log-divergence if we take into accout the weight properly. 
The final log-divergence can be renormalized into the boundary parameters I (or T) and u. See 
SecEJ 

6. Meaning of Introduction of the Weight Function (1): Minimal Area Principle 

In Ref.|19j. we presented the following way to determine the form of W. Let us explain it in the 
warped case. W(p, z) appears as 



2vr 2 f , fVT 



(2vr)4 A/, 



dp dz p s W{p,z)F T (p,z) . (16) 



(2^) 4 y *ji /u 

The integral region is the rectangular of FigllOl We express the above expression by the following 
path-integral. 



r r l / T 
E% u (u,,T)= Vp{z) / dzS[p( 

J Jl/u 



2vr 2 

S\p(z),z] = j— i p(zfW{p(z),z)F^(p(z),z). (17) 



All possible paths {p(z)\uj 1 < z < T 1 } are summed. The dominant contribution, in the above 
path-integral, is given by 5S = 0. 

d - d\n(WF) 

Dominant Path p w {z) : ^- = d ? n \ WF) ■ ( 18 ) 

p dp 

This result says the dominant path pw{z) is determined by W(p, z). A concrete example is the 
valley bottom line in Figll2l 

On the other hand, there exists another path which is determined independently of the above 
dominant path. That is the minimal area curve r g (z) which is given by 

4 r"r 1 1 

Minimal Surface Curve rJz) : 3 H — r r 5 = , — < z < — , (19) 

z r + 1 ijJ T 

where r' = jf-jj , r" = . This differential equation is obtained by minimizing the area A of the 
hyper-surface ([9]): 5 A = where A is given by 

*Z/^t£'^ diiC ^ 

ds 2 = (<U + , A2 ) — 2~2~ = 9ab(x)dx a dx b , 
A = / ^det^ d 4 x = -^ T \/ r,2 + 1 ( 2 °) 

Hence this path r g (z) is determined by the induced geometry g a b{x). We require [15] the following 
relation in order to define W(p,z). 

Pw(z)=Pg(z) (=-L) ■ (21) 

In the above procedure, we have defined the integral measure d A pEdzW(p, z) by the 5D bulk 
geometry or by the 4D induced geometry. 

7. Meaning of Introduction of the Weight Function (2): Fluctuation of Space-Time 

For the purpose of naturally introducing the idea of the previous sections, we newly define 
Casimir energy in the higher dimension by the following generalized path-integral, (warped 
case) . 



£cas (w, T, A) = J *dp J ri i /u) ) \{Vx a (z) x 
VA =r(l/T) a ' Z 



F(-,z) exp 
r 



= P 
1 [ X I T 1 



2a' Ji/u uj 4 z 4 



+ 1 r"dz 

l/w 



(22) 



where p = AT/ui is the IR-cutoff parameter, and we take the limit AT -1 — > 00 at the final stage. 
l/2a' is the string (surface) tension parameter (a' has the dimension of [Length] 4 ). F(p, z) comes 
from the quantum fluctuation of the bulk matter field. 

The above path-integral expression says the 4D coordinates x a play the role of operators of 
the quantum statistical mechanics [20] . The extra coordinate z is the parameter of the inverse 
temperature. It looks that the space-time coordinates are fluctuating. 

We expect the numerical or analytical evaluation of Casimir energy (|22p directly gives similar 
results explained so far. 




8. Conclusion 

We show, in Figll3( the Planck's radiation distribution in the stereo-graphical way by adding 
the inverse temperature axis. The behavior looks similar to Figll2l except the sign. It says the 
extra coordinate corresponds to the inverse temperature. 

The log-divergences in Casimir energy ()14p and (|15p are familiar in the quantum field theory. 
For the warped case, they are renormalized into the boundary parameter uj. 

= (l - 4cln(£) - 4e'ln(£)J = -a{u r f , 



u r = u 



Wi_ 4c ln(^)-4c'ln(^) . (23) 



Local counterterms are unnecessary. Divergences are directly absorbed into the boundary 
parameter. Note that c and d are pure numbers and represent the interaction between the 
bounaries and the (free) field. Compare this with the ordinary renormalization, like QED, where 
the /3-function depends on the coupling. When c and c' are sufficiently small, the /^-function is 
given as 

|c|<l,|c'|<l , u) r = w(l-cln(A/o;) -c'ln(A/T)) , 

a(mA) lo 

The scaling behavior of u is determined by the sign of c + d . Because we identify Casimir 
energy ([23]) with the cosmological constant, the sign also determines the scaling behavior of the 
constant. Note that the direction of flow, which determines the attractive or repulsive force, is 
not given by the derivative (w.r.t. the boundary parameter) of Casimir energy. 

For the flat case (|14p . the other boundary parameter / is renormalized[15j. 

Parameters, which appears in the 5D warped model, can be fitted in the way consistent with 
the present observation (except the sign). First we take the 4D momentum cut-off A = M p i ~ 



1.2 x 10 19 GeV. The present final result {23]) says |S(Euclidian Action) | ~ / d A XyJ—g\ b s /GN oc 
to 4 . The observational data says ^(Euclidian Action) | ~ / d A x^—g\ b s /GN ~ 
(i? cos ) 4 (10 _3 eV) 4 where R cos ~ 5. x 10 41 GeV _1 is the size of the present universe. The 
experimental result about the Newton's gravitational force tells us the warped parameter oj 
is taken as oj ~ 10~ 3 eV. Note that oj ~ yj M p i / R cos , R cos ■ oj ~ y/ M p i ■ R cos = \/Ndl- Hence 
^(Euclidian Action)] ~ N DL 2 = 4. x 10 121 . From {23]), \E Caa \ ~ AT _1 o; 4 . We can identify 
T- 4 \E Cas \ ~ AT" 5 o; 4 with ^(Euclidian Action)] ~ {R cos ) a oj a . Hence T ~ Rcos^Ndl 1 ^ ~ 
3. x l(T 30 GeV. \i = AT/oj ~ M pl N DL - 3 / 10 ~ 7. GeV. The IR cutoff /i is near to the nucleon 
mass or the weak boson mass. If we interpret, in Fig[8l the number of small 4D balls within 
the S 3 boundary is the degree of freedom of the present system, then it is grossly given by 
A 4 /// = oj a /T a ~ N DL 6 / 5 ~ 9. x 10 72 . Finally we note the neutrino mass, m u , is similar to the 
warp parameter (5D bulk curvature) oj. 

As stated in SecJTl the present formulation of the higher dimensional quantum field theory 
gives the picture of the fluctuating space-time. It is known, in the string theory, the uncertainty 
principle appears in the space-time coordinates [21]. 
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